The T = heavy fermion quantum critical point as an orbital selective Mott transition 
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We describe the T — quantum phase transition in heavy fermion systems as an orbital selective 
Mott transition (OSMT) using a cluster extension of dynamical mean field theory. This transition is 
characterized by the emergence of a new intermediate energy scale corresponding to the opening of a 
pseudogap and the vanishing of the low-energy hybridization between light and heavy electrons. We 
identify the fingerprint of Mott physics in heavy electron systems with the appearance of surfaces 
in momentum space where the self-energy diverges and we derive experimental consequences of this 
scenario for photoemission, compressibility, optical conductivity, susceptibility and specific heat. 
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Heavy fermion materials containing electrons in open 
4/ or 5/ shells and in broad spd bands, continue to be a 
subject of great interest in condensed matter physics 
The description of the antiferromagnetic-paramagnetic 
transition in these systems is highly non trivial, because, 
in addition to the fluctuations of the magnetic order pa- 
rameter, one has to take into account the changing char- 
acter, from itinerant to localized, of the / electrons^, Q. 
Recent publications 0,0] have debated whether the quan- 
tum phase transitions observed in heavy fermions can be 
described as an orbital selective Mott transition (OSMT), 
i.e. a Mott transition taking place in the / orbitals with 
the spd orbitals remaining itinerant. Some finite tem- 
perature aspects of this phenomena are captured by sin- 
gle site dynamical mean field theory (DMFT)^, i, 0], 
a method that captures well the peculiar dynamics of 
the Mott transition assuming a purely local approxima- 
tion. In Refs.[^,[l,0| however, the authors have shown by 
means of slave boson techniques that going beyond the 
local approximation is an essential ingredient to obtain 
an OSMT at zero temperature. 

In this paper we overcome the spatial limitation of 
DMFT by using one of its cluster extensions (the cellular 
DMFT, CDMFT[i3]), which takes into account short- 
ranged correlation, and we demonstrate the existence of a 
T = OSMT. In Ref.[lH, we have fully characterized the 
phase diagram of a heavy fermion model across a quan- 
tum critical phase transition, separating a strongly renor- 
malized Fermi liquid from an antiferromagnetic phase. 
Here, by constraining the mean-field non-ordered solu- 
tion, we focus on the qualitative evolution of the elec- 
tronic structure. In this way, we isolate the physics that 
stems directly from the localization of the / electrons 
from the physics of the magnetic order that intervenes at 
low temperature in a given material. 

We show in particular that at the transition a new en- 
ergy scale emerges. In this energy range a pseudogap 
opens in the / spectra and the hybridization between 
heavy / and light spd electrons goes to zero, leading to a 
complete decoupling of the two bands. Beyond this en- 



ergy range the f-spd hybridization remains finite. These 
phenomena have a clear interpretation in terms of an 
OSMT, revealed by the appearance of surfaces of diverg- 
ing self-energy in momentum space, fingerprint of a Mott 
mechanism (mottness[l3])- In the conclusions, we derive 
a set of experimental consequences relevant for the nor- 
mal state of real materials close to the quantum critical 
point at temperatures above the ordered state. 

We study the quantum phase transition driven by a hy- 
bridizing parameter V in the periodic Anderson model, 
which describes free spd electrons locally hybridized to 
non dispersing strongly correlated / electrons. The 
Hamiltonian is: 

H=J2ieu- ^^)dlJ^^ + VY, [flj^^ + h.c 

k k 

k i 

where d\^^ [/^^] creates an spd [/] electron with momen- 
tum k and spin a. The conduction band dispersion is 
Ek = — l/3(cosfc:r + cosfcj^ -I- cosfcz), the other parame- 
ters J7 = 10, /i = 0.2 and Ej - ^ = -5.7. The spd and 
/ electrons Green's functions can be written in terms of 
the / electron self-energy E: 



(2) 



where a = f,spd, Xf{k,uj) ^ fi — Ef — I](k, w), 
Yf{k, w) = ^ - £k and X^pd = Yf, Y^pd = Xf. 

We implement CDMFT on a two-site cluster [lH. We 
believe this is the minimal unit able to capture the 
physics close to the transition point. The Hamiltonian 
in Eq. ([1]) is mapped onto an effective two impurity An- 
derson model (2IAM) and solved self-consistently via the 
Lanczos method [isj. which introduces a finite energy res- 
olution on the Matsubara axis[l3| w„ = (2n— l)7r//3, with 
(3 = 100. 

In order to physically interpret our results, we extract 
the momentum dependent lattice self-energy E(k, w) in 
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FIG. 1: The local density of states (DOS) as a function of the 
hybridization V m a. low energy window around the chemical 
potential (Hubbard bands are out of the picture at a; ~ ±5). 
Black continuous lines are the cluster /-electron DOS. The red 
dashed [green dot-dashed] lines are the /[spd] -electron DOS 
~h X/k ^f[apd]i}^,^)- The bottom panels show the effective 
hybridization ReG/_3pd('i', k) (blue continuous line). 



Eq. ^ from the cluster quantities, restoring the cu- 
bic lattice symmetry. Various methods have been pro- 
posed in the literature[14]. Close to the transition, 
where particles tend to localize, it has been shown^lB'l 

that a suitable quantity to adopt is the cluster cumu- 

I -1 

. In our case 



Ef) 1 - 



I- i Mi(cosfc^ , 
M22 and Ml : 
- Ef - M(k,w) 



cos ky + COS kz) 
-- M12 = M21, 
A stringent 



lant M{uj) = {lo + fi 

we have M(k) = Mq 
where Mq — Mn = 
and E(k, w) = uj + ^ 
self-consistent test of this periodization can be obtained 
by re- calculating the local / electron Green's function 
J2kC! f{i<.,uj), via Eq. ([2|), and confronting it with the 
cluster counterpart, direct output of the CDMFT calcula- 
tion. In Fig. [T]we show the low-energy imaginary parts of 
the local / Green's functions (the density of states DOS). 
The good agreement between the periodized / DOS and 
the cluster / DOS validates our procedure. Moreover 
we show the DOS for the spd electrons and the effective 
hybridization ReG/^spd — J2i!.^^^f-spdi^^,^)^- These 
quantities demonstrate that, as a function of the tuning 
parameter V, the system undergoes a phase transition. 
The numerical uncertainties become greater near the 
transition, hence we cannot determine whether the tran- 
sition is second order, as predicted in a scaling theory [l^ 
(in which case the best fit of the pseudogap scaling with 
a power law gives an exponent zv ~ 0.33), or first or- 
der as found in a recent Guzwiller treatment Q. For 
V > V* ^ 0.58 the system is in the heavy- fermion phase 
where the / electrons present a Kondo peak at the Fermi 
level w = and take active part in the conduction. The 
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FIG. 2: Evolution across the transition point {V* ~ 0.58) 
of the /- (top row) and spd- (bottom row) electron spectral 
functions along the path X=(7r,0,7r) — > F = (0, 0, 0) — > 11 = 
(tt, TT, tt) in momentum space. The color scale (bottom legend) 
is h = 1.0[2.0] for /[spd]-electrons. 



strong hybridization with the spd electrons is evident in 
the suppression of the spd DOS and in the non-zero value 
of the effective hybridization / — spd close to a; = 0. The 
intensity of the / peak reduces approaching V*, while 
at the same time the spd spectral weight enhances. For 
V < V* the system is in an orbital selective Mott state 
where the / electron spectrum has a gap. The / electron 
spectral weight is not completely transferred from low 
energy to the Hubbard bands (placed around oj ~ ±5) 
but rather to a new intermediate energy scale, giving rise 
to a pseudogap. Within this pseudogap the spd electrons 
recover the free band DOS and the effective hybridiza- 
tion is zero. This shows that the spd band at low energy 
is completely decoupled from the / band, but the effec- 
tive hybridization remains active at a finite intermediate 
energy scale. 

We can now display the quasiparticle bands along some 
specific cuts in the w — k space. From Eq. ([2]), we notice 
that Gf transforms into Gspd upon the exchange of Ef — 
E(k,a;) with eu. The poles of the / and spd Green's 
function are therefore the same, provided the ImSk is 
small as in Fermi liquid theory, and the same resolving 
equation is obtained for either a — f, spd: 

uj + fi-Ef-ReT.{k,uj)^ V^/iuj + fi-ew) (3) 

The spectral-weight contribution to the electronic bands 
coming from / and spd electrons are however very dif- 
ferent. In Fig. [5] we show the / (top row) and spd (bot- 
tom row) spectral functions —^IrnGf^spd] (k, w) along the 
path X = (tt, 0, tt) -> r = (o", 0, 0) ^ n = (tt, tt, tt) of 
momentum space, for varying hybridization-parameter V 
(from left to right). For V = 0.67 > V* the band cross- 
ing the Fermi level has predominantly /-character at low 
energy and a strongly renormalized effective mass. In 
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FIG. 3: Evolution of the Fermi Surface (FS, red) and the 
surface of diverging self-energy (ZS, blue) across the transition 
point V . An important role is played by the free conduction 
electron Fermi Surface (FSO, green). The Mott character of 
the transition is marked by the divergence of the self-energy, 
detected by the cluster self-energy eigenvalues T,ev at 1/ ~ 
0.58 and T,c,dd at V ~ 0.43, as shown in the V — i(jj„ space 
(bottom panel). 



approaching the transition point the / electron contribu- 
tion quickly reduces until disappearing completely from 
the Fermi level for V ^ 0.58 = V*. In addition, beyond 
the transition the / band shifts to negative energies. In 
describing the localization of the / electrons therefore, 
the double effect of suppression and translation of the / 
band has to be taken into account. Our result is a predic- 
tion that can be observed in photo-emission experiments. 
At the same time, in crossing the transition, the effective 
mass of the spd electrons is reduced to the free value. 

Recent studies [l5l| have shovifn that insights into quan- 
tum phase transition phenomena can be attained by 
studying not only the Fermi surface (i.e. poles of the 
Green's function), but also surfaces of zeroes in the 
Green's functions (i.e. poles of the self-energy). We first 
remark that in our model there is always a G/ = surface 
in momentum space, corresponding to the free conduc- 
tion electron Fermi surface FSO (given by Ek — /i = 0). 
Further surfaces of zeroes ZS in G/ can appear in k-space 
if there are k points for which E(k, 0) oo. In this case 
we observe that Ggpd reduces to the free conduction elec- 
tron Green's function (Eq. ^ . We show that this latter 
phenomenon indeed takes place in approaching the tran- 
sition point V* . In Fig. [3] we present the Fermi Surface 
FS (determined by w ^ in Eq. O, FSO and ZS for 



different values of the hybridizing coupling V across the 
transition point. For convenience sake, only the lower 
half of the 3-dimensional Brillouin zone is shown. In the 
heavy-fermion phase V = 0.67 > V*, only FS and FSO 
are visible at w = and far apart in momentum space. 
In this case FSO is not relevant for the low-energy physics 
of the system. As soon as V ^ V* , however, a small ZS 
appears around the point k = (0, 0, 0) which pushes FS 
to collapse onto the free FSO. Since at FSO G/ 0, 
this effect originates the strong suppression and disap- 
pearance of the / spectral weight at the Fermi level (see 
Fig. [2]). The appearance of ZS can be already seen in 
the cluster quantities, which are displayed in the V — iu!„ 
space in the bottom of Fig. [3l At V = V* a diver- 
gence takes place for ujn 0[11] in the even eigenvalue 
of the cluster self-energy Set, (left panel), which, via the 
periodization procedure, corresponds to the lattice self- 
energy at k = (0, 0, 0). By further reducing V below V* , 
ZS travels from k = (0, 0, 0) to k = (tt, tt, — tt), where the 
divergence appears in T,odd for V — 0.43 (right panel). 
The position in k-space of the FS remains unchanged 
for V < 0.58, numerically overlapping with FSO. This 
indicates that at the Fermi level G/ = 0, i.e. the / elec- 
trons remain in a Mott state, and Gspd reduces to the 
free Green's function. 

The appearance of a divergent self-energy proves that 
mottness is the physical mechanism governing the lo- 
calization of /-electrons. In an OSMT not all orbitals 
undergo a localization. In the metallic phase all the 
orbitals participate in determining the Fermi volume, 
but, after the transition took place, some "selected" or- 
bitals do not contribute to the Luttinger counting any- 
more. Across this transition a change in the compress- 
ibility of the system is expected as localized orbitals be- 
come incompressible. This is observed in the actinide 
series, where the Mott transition can be driven e.g. by 
pressure A Mott transition is also characterized by 
a significant rearrangement of the electronic structure, 
since there is a transfer of spectral weight from low to 
high energies. In our case the spectral weight is not en- 
tirely transferred from the Fermi level to the Hubbard 
bands, but to an intermediate energy scale giving rise to 
a pseudogap. The resulting modifications of the quasi- 
particle dispersion can be understood in terms of diver- 
gence of the self-energy similarly e.g. to the pseudogap of 
cuprates (see ref. 15|)- Hence, in general, both a Fermi 
volume change and a significant rearrangement of the 
bands are expected when a material undergoes a Mott 
transition. Experimentally this would be detected by 
jumps in the Hall coefficient [3] and in the de Haas- van 
Alphen frequencies [l^. Measuring the phonon disper- 
sions as a function of temperature is another powerful 
probe of the orbitally selective Mott transition, as sug- 
gested in Ref. [13] ■ Furthermore the analysis of the 2IAM 
underlying our self-consistent solution llj suggests other 
experimental predictions. Close to the transition the 
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particle-hole symmetry breaking in the lattice model gen- 
erates a leading irrelevant operator which is forbidden in 
the symmetric case and that causes a log T divergence in 
the specific heat coefficient at high temperature [21[ , to- 
gether with a log T divergence in the staggered spin sus- 
ceptibility and in the pairing susceptibility. On the other 
hand, the formation of singlet correlations on the energy 
scale of the pseudogap in the Mott insulating phase im- 
plies a depression of the uniform spin susceptibility at 
low temperature which has been observed for example in 
CeRhln^ but not in CeC'oIn^ 22, 2^. Since these effects 
originate from the competition of Kondo screening and 
RKKY interaction, we expect them to take place below a 
temperature comparable to (the largest of) these two en- 
ergy scales (roughly 10"'^ the bandwidth, corresponding 
approximately to 10 K in CeRh[Co]In5). This can be 
understood considering that CeRhln^, unlike CeCoIn^, 
is antiferromagnetic at zero temperature and hence lies 
on the Mott selective side of the transition. Optical con- 
ductivity (not shown) displays a clear hybridization gap 
in the delocalized phase which is absent in the phase 
where the / electrons are localized. This is consistent 
with the experimental assignment of CeRhln^ to the lo- 
calized side and CeCoIn^ to the itinerant side of the 
transition 
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To summarize, we have discussed the evolution of the 
momentum resolved spectra of the periodic Anderson 
model across the quantum phase transition, showing its 
orbital selective Mott character. We have described how 
the electronic structure undergoes dramatic reconstruc- 
tion: at the transition the / spectral weight is completely 
suppressed at the chemical potential and a new energy 
scale emerges in the form of an /-electron pseudogap. 
Within this latter energy range, the spd band reduces to 
the free band. / and spd electrons become totally decou- 
pled at low-energy while retaining a finite hybridization 
at higher energies. We have finally shown that the con- 
cept of a surface of diverging self-energy is useful for the 
understanding of this phenomenon. 

It is important to stress that the orbital selective Mott 
phase studied here is not a stable phase at T = be- 
cause the / electrons order magnetically as soon as they 
decouple at low energy from the conduction band[l3|. In 
our calculations the magnetic ordering originates from 
an instability of the orbital selective Mott state and not 
from an instability of the itinerant paramagnetic heavy- 
fermion state. This supports the interpretation of the 
magnetic transition as a byproduct of the OSMT (see also 
Ref.;d|). Our CDMFT study improves previous DMFT 
studies where the local character of the theory for- 
bids the T = OSMT, and previous slave boson studies 
@) B m, Ell > where a finite bandwidth in the / electrons 
must be introduced in order to have an exchange mech- 
anism that is not killed by the vanishing of the effective 
hybridization. In our case such a term is not needed, 
because retaining the full frequency dependence of the 



self-energy allows to have a vanishing effective hybridiza- 
tion at the Fermi level, but at the same time an exchange 
mechanism generated by the non- vanishing hybridization 
at finite frequency. 
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